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ABSTRACT 

Present methods for determining sound speed corrections 
for echo soundings in continental shelf areas are time-consuming 
and expensive. This study was undertaken to determine whether 
or not sound speed correctors of sufficient accuracy could be 
deduced from historical data. 

Historical sound speed data for an east coast, shelf area 
indicates that temporal and spatial variability exceeds 
acceptable limits for sounding corrections, thus precluding 
the use of historical data only for corrector determination. 

Examination of temperature and salinity data indicates 
that historical salinity values in the region are sufficiently 
stable to allow acceptable sound speed corrections to be made 
derived on the basis of in-situ temperature measurement and 
historically derived salinity values. Expendable bathythermo- 
graph probes are capable of meeting temperature requirements. 

The constancy of the T-S relation in such regions may 


allow salinity determination from water temperatures alone. 


z TABLE OF CONTENTS 


- I. INTRODUCTION ------------------ 9 
OA CORNET er pe rs ae ne eee ge 9 
B. ACCURACY REQUIREMENTS - ----------- 11 
C. HISTORICAL TECHNIQUES - - ------+----- 12 
D. MEASUREMENT AND COMPUTATIONAL _ | 
REQUIREMENTS-SOUND SPEED EQUATIONS - - - - - 15 
II. PROCEDURE AND METHODOLOGY = - - = = = su =< 20 
Bh, SRCRY ARES 2 eae See Se ee Re 20 
B. DATA SOURCE - = = + ----+-+--<+-++-+-+--s 30 


C. TECHNIQUES FOR SUMMARIZING 


SOUND SPERD DATA =<. = «45 ee ee uu s 31 

ty. COMPUTER BROCEEaING - = = = = 2 oe we wane 34 

E. BROCRDURE = 2 <2 ce ee cect aweens 40 

Tit. DISCUSSION OF RESULTS = 2 = = ss sew sane 46 
TV. CONCLUSIONS << 4 — ws ue ee meee wen eo 54 
APPENDIX 1: COMPUTER PROGRAMS AND DOCUMENTATION - - - 58 


APPENDIX 2: VELDAT, SALDAT AND TEMDAT RESULTS 


AND SEASONAL PROFILES - - - - - - - - - - 104 
APPENDIX 3: SEASONAL STATION PLOT - - - - - - - - - - 122 
BISLIOGHAPEY «<< Ge eee tne we eee 126 
INITIAL DISTRIBUTION LIST - - +--+ -<«s 2 ssc 129 
FORN DD 1473 «ss eee nen ewe nnens he a 1 


Iv. 


LIST OF TABLES 


Temperature (°C) in the Western North Atlantic - ~ - 
Salinity (°/..) in the Western North Atlantic - - - 
Comparison of HIDAT Results with 

TEMDAT and SALDAT Results - ------7------>- 
Summary of VELDAT RESULTS 

Region l ----+---+-+--e-+7-7+47727-47-72--°7 
Variability of Sound Speed Temperature 


and Salinity -----+-+---s+--e-++7277724-°7 


27 
27 


45 


47 


51 


LIST OF FIGURES 


Schematic Representation of BS? Sounding 
System Showing Swath Coverage 


Index Map of Study Area 

Surface Circulation in Study Area 
R/V GILL Station Plan 

Water Mass Distribution 

Sample NODCRD Output 

Sample Input Data Cards 


Regional Division of Study Area 


ACKNOWLEDGEMENTS 

I would like to express my sincere appreciation to | 
LCDR A. B. Chace, USN, and CDR D. E. Nortrup, NOAA, as thesis 
advisor and reader, for encouragement and guidance throughout 


this project. J 


I am especially indebted to LCDR Alan Pickrell, NOAA, for 
| his invaluable assistance in developing programming for this 
analysis. 

I would also like to thank Dr. R. G. Paquette for his 
kind assistance and permission to utilize existing computer 
programs written by him. Further, I would like to extend my 
appreciation to Mr. Henry Odum of the National Oceanographic 


Data Center for his timely assistance in providing the data 


Finally, I would like to thank my wife for her patience 


f 

tapes used in this study. 

and understanding during this undertaking. 
f 


f 
¥ 
; 
F 
{ 
i 
i 
f 
i 
t 
qj 


I. INTRODUCTION 

A. GENERAL 

The objective of this thesis was the application of basic 
oceanographic principles to the solution of an operational 
hydrographic problem. Depths determined by echo sounder are 
based on an assumed speed of sound in seawater. Corrections 
must be applied to such depths to account for the difference 
between the assumed and actual speed of sound. Present prac- 
tice dictates that in-situ measurement of the sound speed 
profile in a working area be made and depth corrections com- 
puted. The problem is to determine if temporal and spatial 
sound speed variability in continental shelf areas is small 
enough to allow use of historical salinity and temperature 
data archived in the National Oceanographic Data Center's 
(NODC) Oceanographic Station File to correct echo soundings. 
Historical data was examined to determine whether or not 
correctors of sufficient accuracy could be obtained from 
historical data without resort to in-situ measurements. 

Applicability of historical data (in the form of tempera- 
ture and salinity measurements) has been widely discussed but 
limits of applicability have not been well-documented (Mobley, 
1977). 

This study was undertaken to examine the historical data 
available, its variability within a sample coastal region, and 
limitations on application in meeting National Ocean Survey 


(NOS) hydrographic requirements. 


Echo Sounders (or fathometers) do not measure depth 
directly but rather measure the time delay between an outgoing 
sound pulse and the return echo of this pulse. Depth is then 
derived by dividing the round crip travel time by two and 
multiplying this value by the assumed speed of sound in sea- 
water. The fathometer makes this transformation electronically 
or mechanically within the device itself and displays depth 
of water beneath the transducer. The depth displayed is the 
nominal or fathometer depth. True depth can be determined only 
if the soundspeed profile through the water column is known 


and used to calculate depth as in the following equation: 
2At 
fi oa s/f V(Z, )dt (1) 
/0 


"Where V(Z,) is the sound speed at the level Z. where the 
signal passes at the time 0<t<2At. Corrected depth is derived 
from nominal depth by applying a correction according to some 
standardized procedure" (Greenberg and Sweers, 1972). 

Most echo sounders used in hydrographic surveying are 
constructed with a calibrated speed of sound assumed in the 
device. This speed is usually either 1463 meters per second 
(800 fathoms per sec) or 1500 meters per second (820 fathoms 
per second). Reference speeds are chosen to represent 
realistic ocean values which are convenient for design. Echo 
sounders in use by the National Ocean Survey at present are 


calibrated with an assumed sound speed of 1463 meters per 
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second. This value is reasonably near the actual speeds en-~ 
countered in most waters surveyed by NOS (Umbach, 1976). 

Sound speed corrections are applied to bring standardized 
values close to actual vertical column sound speeds for a 
location (Umbach, 1976). 

Field measurements are made during the course of a survey 
in order to determine representative vertical sound speed pro- 
files. These field measurements and their reduction are often 
time consuming, expensive and result in considerable conjecure 


regarding areas and times of corrector applicability. 


B. ACCURACY REQUIREMENTS 


The requirements for accuracy of echo soundings for hydro- 
graphic surveys are based on standards agreed to by member 
states of the International Hydrographic Bureau (IHB). Depth 
measurement requirements are stated in Special Publication 44, 
"Accuracy Standards Recommended for Hydrographic Surveys", 
published by the IHB in 1968. These standards vary with depth 


and may be summarized as follows (Umbach, 1976): 


(a) O-11 fathoms - - ---+--+-- - LiO* cht 
(b) 11-55 fathoms ~ - - - - - - - 3.0 £t. 
(c) Deeper than 55 fathoms - - - - 1% of depth 


These standards represent the maximum allowable errors in 
measurement of depth due to all sources, including variation 
of sound speed in the water column. ‘"Accuracies attained for 
hydrographic surveys conducted by the National Ocean Survey 


shall equal or exceed those given by the IHB (Umbach, 1976)." 


The National Ocean Survey's requirements for accuracy in 
determining the speed of sound for correcting echu soundings 
are stated as follows: ''The velocity of sound must be known 
with sufficient accuracy to ensure that no sounding will be 
in error by as much as 0.25% of the depth from this cause 
alone. Therefore the mean velocity of sound must be known to 
within +4 meters per second" (Umbach, 1976). These standards 
apply only to conventional echo sounding systems and survey 
operations involving these systems. An implicit assumption 
is that the sound pulse or "beam" is vertical. 

Future systems under test by NOS such as the Bs? 
(Bathymetric Swath Sonar System) will require greater accu- 
racies. This is due to the fact that many beams are projected 
at varying angles through the water column to attain swath 
coverage. Figure 1 is a diagram of such a system. The longer 
acoustic path lengths, as well as varying trajectories, makes 
this technique more sensitive to sound speed error than 
present vertical depth measurements. Accuracy requirements 
for such a system have been stated as +2 meters per second 


(Mobley, 1977). 


C. HISTORICAL TECHNIQUES 


A brief review of historical techniques utilized by the 
hydrographer for determination of sound speed is necessary in 


order to properly appreciate the magnitude of the problem. 
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Figure l 
Schematic Representation of BS’ 
Sounding System Showing Swath Coverage 
(From Mobley, 1978) 
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There have been at least four common methods used to 
determine either the sound speed profile or the correction in 
depth units to be applied to soundings (Maunsell, 1976). 

1) Bar checks or leadline comparisons in which the depth 

recorded by the echo sounder is compared to the depth 
of the reflecting surface of the "bar" lowered to a 
known depth, or the leadline depth to the bottom. 
This procedure has usually been limited to use in 
shallow depths due to problems of manuevering the 
bar under the transducer, and correct gauging of the 
depths indicated by the leadline. This method gives 
a correction to be applied to the nominal depths in- 
dicated by the echo sounder. Sound speed may be 
determined from depth differences provided that 
draft of the transducer is known. 

2) Direct measurement of salinity temperature and depth 

by either Nansen casts, salinity, temperature and 


depth sensors (STD) or the expendable STD (XSTD). 


Sound speeds at various standard depths are then 


computed using an equation relating salinity, tempera- 


ture and pressure to speed, such as Wilson's 1960 
equation. 

3) Direct measurement of sound velocity as a function of 
depth by lowering a sound velocimetwr to various 
depths and recording the speeds obtained at these 


depths. These devices have been poorly utilized in 
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the past due to design limitations. The development 


of an expendable Sound Velocimeter (XSV) makes this 
device promising for the future. 

4) Historical tables or atlases of sound speeds or salinity 
and temperature data observed in the same area and 
season as the operations. Tables such as Matthews 
Tables (Matthews, 1939) and Heck and Service's Tables 
(Heck & Service, 1924) have been in use for quite some 
time with periodic efforts being made to improve upon 
the accuracy of results obtained (Ryan, 1974). These 
tables have not been generally applied to depths of 
200 meters or less (Sherwood, 1974). 

In general then, there have beena variety of methods employ- 
ed, some utilizing direct measurement and others utilizing 
historical information. All of the methods summarized have 
been noted for deficiencies or disadvantages of one form or 
another. Adequate direct measurements of the sound speed pro- 
file in the operational area provide the best accuracy; 
however, in terms of expense and efficiency, there may be 


useable alternatives. 


D. MEASUREMENT AND COMPUTATIONAL REQUIREMENTS - SOUND SPEED 
EQUATIONS 


Sound speeds in a fluid may be calculated with various 
equations relating speed to density and elasticity of the 


medium. The following equation extracted from ''Fundamentals 


of Acoustics" by Kinsler and Frey relates speed of sound in 


a fluid to the isothermal bulk modulus (Br), the ratio of 
specific heats (I) and density in kg/m° (P4)- Isothermal 
bulk modulus is an elastic modulus measuring the difficulty 


of compressing a liquid (Kinsler and Frey, 1962). 


c ay! Bp (2) 


Po 


In the case of seawater, these properties are functions 
of Salinity, temperature and pressure (Ryan, 1974). 

Standard practice within the National Ocean Survey has 
been to use a form of Wilson's 1960 equation for the speed 
of sound as a function of salinity, temperature and pressure 
(Umbach, 1976). The usual case has been to calculate the in- 
situ speed profile from oceanographic observations within a 
project area. This empirical equation was developed as a 
"best fit'' to certain measured data. The 1960 equation was 
an improvement over earlier versions in order to expand the 
range of applicability. The equation is stated by Wilson 
to fit data within the ranges of temperature from ~4°c to 
30°c, salinities from 0 °/.., to 37 °/,. and pressures be- 
tween 1 kg/cm and 1000 kg/cm?. This equation is said to 
have achieved a standard deviation of 0.30 meters/second 
from the mean for all data obtained (Wilson, 1960). 

It has been noted by other researchers that significant 
accuracy improvements have been achieved in determining sound 


speeds utilizing equations that fit data limited to real 
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seawater conditions (LeRoy, 1969). One of the two equations 
formulated by LeRoy in 1968 "fit Wilson's second equation for 
seawater to within 0.1 meter/second in the domain described 
and the other fits Wilson's corresponding data with a better 
accuracy than does Wilson's equation" (LeRoy, 1969). Other 
equations such as those formulated by Frye and Pugh in 1971 
also claim similar improvement (.lm/sec error estimate) 

(Frye and Pugh, 1971). 

These variations are discussed as a means of pointing 
out possible error sources and improvements to techniques to 
achieve desired accuracies in sound speed determinations. 

Wilson's 1960 equation is stated in the Journal of the 


Acoustical Society of America, October 1960 as follows: 


V=1449. 14+V,+Vp+V5+Vonp 

Vp24.5721T-4.4532X10~ “17-2. 6045x10-*T? 
+7.9851x10~°T# 

Vp=1.60272K10~ "P+1.0268K10~ °P” 


-9,3 ~12)4 


+3.5216X10 ~P°-3.3603X10 (3) 


Vg=1.39799(S-35)+1.69202K10" 9(s-35)? 


2 2 


Voqp* (8-35) (-1.1244X10" T+7.7711X107 'T 


+7.7016X10~ °p-1. 2943x10 ’Pp?+3. 1580x107 Spr 


~9pq2 )+p(-1. 8607X107 47 


+7. 4812x107 °T+4 5283x107 81?) 


7r+1.8563X10~°T7) 


+1.5790X10 


+P?(-2.5294x107 


+p?(-1.9646x107)°T). 


In this equation the units of temperature, pressure, 
salinity, and sound speed are °C, kg/cm, parts per thousand, 
and meters/sec, respectively. This equation is essentially 
the same as that given in the NOS Hydrographic Manual with 
the exception of the term containing variation of gravity 
with latitude. This term is used if pressure is not measured 
but computed from depth. The equation as stated in the Hydro- 


graphic Manual is: 


Vy = 1449. L4tVptV tVotVntVorp (4) 


"Where *e is the speed of sound in meters per second, Vp is 
a correction for pressure, Vd? is a correction for variation of 
is 


gravity with latitude, Vg is a correction for salinity, Vip 


a correction for temperature and V is a correction for the 


STP 
combined effect of salinity, temperature and pressure” (Umbach, 
1976). All terms in equation 4 are the same as those in 
equation 3 with the exception of the "s term. Ya may be deter- 
mined by reference to tabulated values in Table 12C of Special 
Publication 68,''Handbook of Oceanographic Tables" published 

by the U.S. Naval Oceanographic Office in 1966. 

It has been necessary to examine these equations in order 
to determine the effects which small variations in temperature 
or salinity have on the results attained. 

As stated in the Hydrographic Manual, utilizing Wilson's 


Equation, temperature measurement accuracies of +1°C and salinity 


measurement accuracies of +lppt are required to satisfy sound 


speed correction requirements (Umbach, 1976). 
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An analysis of the Wilson equation was performed by 
Testing Division, Office of Marine Technology (OMT), NOAA in 
1976 to determine sensitivity of the equation to variations of 
temperature and salinity (Bivins, 1976). Their results were 
stated as follows: "if corrections were made on the basis 
of temperature measurement with a one sigma (lo) accuracy of 
+.1°c, a natural variability of +3ppt (lo) from the measured 
salinity could be tolerated and the NOS sound speed accuracy 
of +4 meters per second would be satisfied" (Bivins, 1976). 
This analysis pointed out that the most critical measurement 
was temperature and if temperature accuracies of .1°C were 
met then salinity measurement errors within a larger range 
than previously noted could be tolerated. 
The accuracy requirements, as well as the OMT analysis 
all bear on the acceptability of historical information for 
this application. In fact, the requirements stated are the 
criteria for acceptance or rejection of the technique proposed. 
Sound speed corrections are computed within NOS by the 
"Summation of layers'' method (Umbach, 1976). Once temperatures 
and salinities have been determined at standard depths, sound 
speeds at mid-layer depths are computed by graphical or 
numerical methods. These layer corrections are then summed to 
produce the total sound speed corrections applicable to given 
depths (Umbach, 1976). Layer thicknesses utilized in this 
procedure are specified in the "NOS Hydrographic Manual" as 


10 meters for depths to 200 meters. 
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II. 


PROCEDURE AND METHODOLOGY 


A. STUDY AREA 

The area selected for study was a coastal region off-shore 
of Charleston, South Carolina. It was originally decided to 
examine a one-degree "Square'' extending from latitude 32°N 
to latitude 33°N and longitude 79°W to longitude 80°W. The 
study area was subsequently enlarged slightly in order to add 
additional data points near the northeast corner of the one- 
degree square. Figure 2 is an index map showing the study 
area. 

The area was chosen on the basis of several factors: 

First, it appeared to be representative of the oceano- 
graphic environments characteristic of east coast survey areas. 
The influx of estuarine waters in the nearshore area, bottom 
topography, depth range and other physical features are similar 
to those noted for shelf waters from Georgia to the Virginia 
Capes. 

Second, the size of the region selected was typical of 
a one-year hydrographic project by a single coastal survey ship, 
or a six-month project by two ships. 

Third, examination of the Environmental Data Service's 
"Key to Oceanographic Records Documentation No. 2, Temperature, 
Salinity, Oxygen and Phosphate in Waters off United States" 
showed the data density in this region to be representative 


of that in nearshore regions along the East Coast. Some 300 
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Nansen cast stations dating from 1966 to 1975 were included 


within this area. 

Finally, the area also corresponded to NODC data format 
and limitations of existing computer programs for abstracting 
and manipulating this data. 

The area extends from shore approximately fifty-five NM 
offshore (to the 200 meter curve). The majority of the region 
has a depth range of eighteen to seventy-three meters (10 to 
40 fathoms). Bottom topography is gently sloping until the 
shelf break is encountered. 

Surface circulation is seasonal in character and is in- 
fluenced by the following factors: river run-off, horizontal 
temperature gradients, wind and coriolis force effects on 
tidal motions in restricted waters (Bumpus and Lauzier, 1965). 

The region is mainly inshore of the Gulf Stream although 
the Gulf Stream affects its circulation patterns (Gaskell, 
1972). No systematic study of seawater movement has been 
made, though several generalizations have been based on 
limited data from drift bottles (Kuroda and Marland, 1973). 

The surface circulation has been characterized as having 
a prevailing northerly drift during the winter and spring 
(Bumpus and Laugier, 1965). The pattern becomes less well- 
defined during the summer. During the summer transitional 
period there is water movement both to the south and to the 
north. The northerly movement is apparent in the very near 
shore region and over the outer edge of the shelf. Between 


these two motions is a southerly flow (Bumpus and Lauzier, 1965). 


During the autumn, most of the region is under the influence 
of a southerly drift, with the exception of the outer portion 
of the shelf, where the flow is toward the northeast (Bumpus 
and Lauzier, 1965). Current speeds are maximum during the 
winter and autumn and a minimum during the summer months (U.S. 
Naval Oceanographic Office, 1965). 

Kuroda and Marland (1973) noted that the surface circula- 
tion in this region is strongly dependent on the prevailing 
winds. 

Figure 3 is a diagram exhibiting the seasonal nature of 
the surface currents in the area extracted from the American 
Geographical Society's "Serial Atlas of the Marine Environment, 
Folio 7". 

Systematic study regarding water mass distribution or 
temperature and salinity distribution is limited in this region 
(Kuroda and Marland, 1973). The Bureau of Commercial Fisheries 
R/V THEODORE N. GILL studied this region in 1953 (Kuroda and 
Marland, 1973). Kuroda and Marland (1973) in "Physical and 
Chemical Properties of the Coastal Waters of Georgia" summarized 
the GILL data. Figure 4 is a station plan for the GILL cruises. 

This data shows that the isotherms trend northeast and 
southwest, paralleling the coastline and the approximate axis 
of the Gulf Stream. Monthly minimum temperatures were observed 
inshore and maximums generally seen in the offshore areas. 
Minimum surface temperatures (10°-23°C) occurred in January and 
February and maximum temperatures (27°-29°C) occurred in June 


and July (Kuroda and Marland, 1973). Minimum and maximum 
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Figure 4 


R/V GILL Station Plan 
Cross-hatched Region is Area of Present Study 


(From Kuroda and Marland, 1973) 
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temperatures at depth occurred gradually later with a direct 
depth relation (Kuroda and Marland, 1973). As an example, 
minimum temperature at 50 meters may occur in February whereas 
at 100 meters it occurs in March. The temperature at 200 
meters at a point is nearly constant throughout the year. 
Table I, extracted from Kuroda and Marland (1973) summarizes 
the temperature distribution. 

Isohalines are also shown to parallel the coastline (Kuroda 
and Marland, 1973). Near shore, surface salinities exhibited 
a seasonal change whereas offshore salinities showed little 
seasonal effect and remained near 36°/.,, (Kuroda and Marland, 
1973). 

Near the coast,the minimum salinities of 32°/,, were ob- 
served in January and February, and the maximum value of 35°/,., 
occurred in summer and autumn (Kuroda and Marland, 1973). Of 
particular interest are the near homogenous salinity values 
(36°/,..) observed throughout the area at a depth of 50 meters 
(Kuroda and Marland, 1973). 

Ranges for salinities at greater depths (100 to 200 meters) 
were also small (35°/,. to 36°/,.,.) (Kuroda and Marland, 1973). 
Table II, also extracted from Kuroda and Marland, is a summary 
of the salinity variations observed. 

Analysis of the Temperature-Salinity diagrams for the R/V 
GILL data by Kuroda and Marland revealed the presence of four 
different types of water masses in the region. 


The water mass present along the outer edge of the shelf 


is equivalent to the North Atlantic Central Water (NACW) as 


TABLE | 


TEMPFRATURE (°C) IN THE WESTERN NORTH ATLANTIC 


tre OLE PTH | 
QO METERS 10 METERS 50 METERS | 100 METFRS | 206 METERS 


JAN, - FEB, 


FEB, = MARCH 


TABLE It 
SALINITY (°/oo) IN THE WESTERN NORTH ATLANTIC 
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described by Stommel (1965) (Kuroda and Marland, 1973). This 
type exhibited the highest temperature and salinities observed. 
The water temperature above 100-200 meters showed large sea- 
sonal change but relative constancy below this depth (Kuroda 
and Marland, 1973). 

A type of water labelled "mixing water" by Kuroda and 
Marland (1973), exhibited very similar properties to the NACW 
previously described but temperatures were approximately 5°c 
lower in winter and autumn. The salinities observed above 100 
meters had almost identical values as the first type while at 
200 meters there was a decrease of about .4°/,, (Kuroda and 
Marland, 1973). The third water mass, described as shelf 
water, was nearer shore. Salinities between 34.5°/,, and 
36°/.. were characteristic and exhibited much seasonal change 
(Kuroda and Marland, 1973). 

A fourth water mass was identified as "coastal water" and 
occurred very near the shore. This water mass exhibited large 
seasonal changes of both temperature and salinity (Kuroda and 
Marland, 1973). Convenience in analysis of regional conditions 
was the primary purpose of Kuroda and Marland's classification. 
Figure 5 is a generalized map of the distribution of these 
water masses, extracted from Kuroda and Marland (1973). 

The suite of water masses described extend out to approxi- 
mately 100 miles offshore, just inside the western edge of the 
Gulf Stream. Each mass is continually interacting with the 


others resulting in modification to the original character. 
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Figure 5 


Water Mass Distribution 
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Seasonal distributions of these masses and their movement is 


also described by Kuroda and Marland (1973). 


B. DATA SOURCE 


The data utilized in this analysis was obtained from the 
National Oceanographic Data Center (NODC) in October 1978. 

A magnetic tape of all oceanographic station data (Nansen cast 
and STD stations) for Marsden square 116 was furnished. A 
total of approximately 300 stations were within the area of 
interest. The data consisted of temperature, salinity and 
sound speed (computed by Wilson's equation) for each station. 
Additional information included: position, date, time, ship 
identification and other meteorological chemical and optical 
information wherever it was measured. 

The magnetic tape furnished was in standard NODC format: 
nine-track, 800 bits per inch, 80 character, full-blocked 
records. Each station record consisted of a master record 
containing positional and identification data and a detail 
record containing the salinities, temperatures and sound 
speeds obtained for each standard depth (to the maximum depth 
of the cast). Standard depths for oceanographic observations 
are a matter of agreement within the oceanographic community. 
NODC files contain records (for depths to 400 meters) at the 
following standard depths (in meters): 0, 10, 20, 30, 50, 
75, 100, 125, 150, 200, 250, 300, 400 (NODC, 1974). Where 
observations do not conform to the standard depths, inter- 


polations are made utilizing a 3-point Lagrange interpolation 
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equation (NODC, 1974). Wherever salinity, temperature and 
depth data are present and the values are not doubtful, sound 


speeds are computed (NODC, 1974). 


C. TECHNIQUES FOR SUMMARIZING SOUND SPEED DATA 


The objectives were to examine and summarize the historical 
sound speed data available and evaluate its suitability for 
determination of correctors for echo soundings in the study 
area. 

Temperature and salinity variability was also examined 
and its relation to sound speed variability noted. Several 
methods for summarizing the sound speed variability in an 
oceanic region have been devised by other researchers. 

Russell (1975) developed a method for selecting typical 
sound speed profiles from oceanographic station data with 
application to interpretations of acoustic system performance. 
The data was abstracted in groups representing estimates of 
natural temporal and spatial populations. The study utilized 
NODC data as well as available XBT data. In this method, the 
ocean regions represented natural oceanographic domains. Each 
domain exhibited similar physical properties and current 
regimes. 

Russeli's (1975) study is limited in application. The 
primary reason for this is difficulty in accessing data bases 
other than those maintained by NODC. The study is further limited 


because it is directed to deep water horizontal propagation. 


The profiles generated are useful as pilot information but 
the variability of these profiles within a region is not de- 
fined with the precision required for application to the echo 
sounding problem. 

Audet and Vega (1974) developed a similar technique to 
that of Russell (1975). The data utilized was exclusively 
the NODC Oceanographic Station File. 

In this technique, all profiles for a region are examined. 
Three profiles which exhibit the smallest variations relative 
to all other profiles are selected, the profile exhibiting the 
least variability being selected as the model profile. The 
other two profiles are compared to the model and the standard 
deviation computed between the model and these two representa- 
tive profiles. This study is of interest due to the methods 
employed. However, the study was primarily concerned with 
deeper oceanic regions and application to variable continental 
shelf areas is inappropriate. 


The HIDAT program in use by the U.S. Naval Fleet Numerical 


Weather Central is used to characterize sound speed conditions 
in ocean regions for forecasting performance of military acoustic 
systems. 

Output for the HIDAT program is a month-by-month summary 
showing the temperature and salinity profile which is found 
to be representative of a region. The representative profile 
is an actual observed profile which shows the smallest standard 


deviation from all other profiles in the region. 
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Standard deviations at standard depths are computed and 
represent variance of the representative profile from the 
four most similar profiles in the area. 

The smallest region for which HIDAT information can be 
generated is a one-degree latitude-longitude square. 

The data base for HIDAT includes all NODC data as well 
as extensive XBT data from other sources. This results in much 
greater data density for an area. The HIDAT program is 
similar to the system suggested by Audet and Vega (1974). 


All three prior techniques were developed in response to 


specific operational requirements. None of the three was 


designed appropriately for echo sounding correction. Each 
technique was concerned with selection of an actual profile 
as the representative profile in a region. This is an impor- 
tant design criteria in that it permits identification of the 
sonic layer depth. Identification of the sonic layer depth 
is important for horizontal acoustic propagation. This was 
the application for each of the three studies examined. 

The selection of an actual observed profile is unnecessary 
for vertical echo sounding because minor chang:s in location 
of the sonic layer depth have little effect on mean vertical 
sound speeds. Further research and adaptation of these pro- 


grams to the echo sounding problem might result in useful 


techniques. 


33 


D. COMPUTER PROCESSING 


All computer processing in this analysis was accomplished 
utilizing the IBM 360/67 computer system of the W. R. Church 
Computer Center at the U.S. Naval Postgraduate School (NPS). 

In addition to the programs written to perform the required 
analysis of NODC data, several existing programs for abstracting 
and plotting NODC data were employed. These programs were 
written by Dr. R. G. Paquette of the Department of Oceanography, 
Naval Postgraduate School. All programs were written in 
Fortran and complete listings and documentation are contained 
in Appendix 1. 

Initially, all station records from the NODC unlabelled 
magnetic tape were transferred to a standard NPS labelled tape. 
This resulted in a faster retrieval time and fewer errors 
in reading the tapes. Although the NODC tapes were compatible 
with tape drives at the NPS computer center, reading provlems 
may occur whenever magnetic tapes are read on a drive different 
from the machine on which they are created. Transfer to 
the NPS standard labelled tape eliminated the errors associated 
with this problem. 

A computer listing of all NODC station header records 
for Marsden Square 116 was made. This listing was produced 
using program NODCDUM, written by Dr. R. G. Paquette of the 
Naval Postgraduate School. These records were examined to 
locate all stations within the study area. The stations of 
interest were all located within two adjacent blocks of 


record numbers. This resulted from the format of NODC data 
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tapes. Within a Marsden Square, further blocking or sub- 
division of records is made by one-degree latitude and longi- 
tude squares. 

A complete listing of the header and detail record for 
each station within the study area was then produced on a 
month-by-month basis. As an example, all stations observed 
in a particular month (all years) were listed together for 
ease in indexing. The program used to make this listing was 
NODCRD, also written by Dr. Paquette. The record produced 
was a hard copy printout of all stations observed each month. 
The information in this listing consisted of both the header 
record (position, date, year, country and ship identification 
and NODC index numbers) and the detail record giving the type 
of observation (Nansen or STD), temperature, salinity, density 


(o,), oxygen content and sound speed at each standard depth 


t 
to the limit of observation. This listing served a reference 
and index function throughout the study. A sample record for 
one station is shown in figure 6. 

Four punched data cards were produced for each station. 
The first card was produced using NODPUNA. It contained 
latitude, longitude, NODC reference number, month in which 
the observation was made, and the sound speeds for each 
depth to the limit of observation or 400 meters, whichever 
was shallower. Speeds were recorded to tenths of meters per 
second. The second card contained the same positional and 


index information but observed temperatures were abstracted 


to two decimal places instead of sound speeds. This card was 
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produced using NODPUNC. The third data card was also similar 
but salinities to two decimal places were abstracted for each 
standard depth. NODPUNB was used to produce these cards. The 
fourth data card was produced as input to a station plotting 
program, CHARTPLT. This card contained latitude and longi- 
tude of the station. The card was produced using NODPUN. 

This additional card with positional information was produced 
in order to avoid programming changes to accomodate a new 
format. Samples of all four cards are shown in figure 7. 

The programs used to produce the first 3 cards (NODPUNA, 
NODPUNB, and NODPUNC) differed only in the record abstracted 
(speed, temperature or salinity). The temperature, salinity 
and sound velocity cards were used as input to the statistical 
program written to summarize the data. 

VELDAT was the basic statistical program written to summa- 
rize the sound speed data. Slight modifications to this 
program, particularly in input format, resulted in programs 
SALDAT and TEMDAT. All three have the same basic functions, 
and produce the same results. VELDAT was used to analyze 
sound speed data, SALDAT for salinity data and TEMDAT for 
temperatures. 

Input to these programs were the appropriate data cards 
previously described and shown in figure 7. 

VELDAT output consisted of the mean sound speed profile 
(at standard depths), the standard deviation at each standard 


depth, number of samples included in the calculation and the 
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location and values of maximum and minimum sound speeds. 
Temperature and salinity data were similarly analyzed. 

The values for sound speeds at each standard depth for 
each station in the study area were input into VELDAT, grouped 
by season. As an example, all stations observed in January, 
February or March are examined together. Other seasons were 
defined as: Spring (April, May, June); Summer (July, August, 
September), and Autumn (October, November, December). This 
seasonal breakdown followed the practice of other researchers 
in the area (Kuroda and Marland, 1973). 

Initially, sound speeds for each standard depth for all 
stations in the area were averaged and the sample mean computed. 
The following equation represents the computation made in 


VELDAT : 


Vv, = 4, (5) 


SV was the mean sound speed at a particular standard depth 
(d) over the entire area. )sVg was the sum of all sound 
speed values at a particular viutiaie depth, and N was the 
number of sound speed values found for a particular standard 
depth. Where no value was found for a particular depth or a 
zero value was encountered, this sample was not included and 
N was not incremented. 


A sample standard deviation was computed for sound speeds 


at each Standard depth utilizing the following equation: 
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S) is the standard deviation of the sound speeds at a partic- 


ular standard depth. SV was the mean sound speed computed 


d 
for the particular depth in question and N was the number of 
samples. The same procedure for dealing with blank or zero 
values was used in this computation as in the computation for 
mean sound speed. TEMDAT and SALDAT utilized the same equations 
for producing mean temperatues and salinities as well as 
standard deviations. Appendix 1 documentation includes input 
or data format required, variable specifications for different 
regions and modes of operation possible. 

The procedure utilized to compute mean sound speeds and 
standard deviations at standard depths was suggested by the 
fact that present NOS practice dictates use of the "summation 
of layers'' method for determination of sound speed corrections 
at depth. This method has been previously described in section 


I.C. of this thesis. 


E. PROCEDURE 


Sound speed data cards for all stations were produced and 
sorted by month. Variability of sound speeds over the entire 
area on a monthly basis was examined initially. Lack of sub- 
stantial numbers of observations during any particular month 


precluded use of this technique for drawing meaningful conclu- 
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sions. 
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Three month seasonal groupings were then combined and 
analyzed using VELDAT. The program output exhibited a larger 
variability than could be tolerated under the accuracy require- 
ments stated. The question of applicability of historical 
data would have been quickly resolved had the variabilities 
over the entire area been small enough to match the stated 
accuracy requirements. However, this was not the case, the 
varibility was too great; thus leading to a form of spatial 
subdivision. 

The locations of stations exhibiting the maximum and 
minimum sound speed values revealed that grouping by area 
might reduce the observed standard deviations. These results, 
examined in view of the variation of physical properties of 
the waters (Kuroda and Marland, 1973), suggested a sub-division 
of the area into two regions. Since sound speed values should 
be similar for waters with similar physical properties, the 
objective of the sub-division was to separate the area into 
two natural populations of sound speed profiles. This method 
followed the procedure cited by Russell (1975), although his 
study was concerned with sub~division on a much larger scale. 

The division was made along a northeast - southwest trending 
line approximately paralleling the 18 meter depth contour. 

The near-shore region resulting from this division was expected 
to exhibit much greater variations in sound speeds and hence 

a greater standard deviation. The offshore region was expected 
to show smaller variations of sound speed values. Figure 8 is 


a map of the area showing the division made. Sound speed 
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analysis was then achieved for each of the regions using VELDAT 


and the seasonal breakdown previously cited. A similar pro- 
cedure was carried out in order to examine temperature and 
salinity variability in the region. TEMDAT and SALDAT were em- 
ployed to this purpose. 

Graphic profiles were produced from the VELDAT output 
showing the mean sound speed profile, and the maximum and mini- 
mum profiles for the offshore region. Profiles were produced 
for each season. Salinity and temperature profiles were also 
produced on a seasonal basis. Profiles for the inshore region 
were not produced due to the sparseness of data in the region. 

Output from the statistical programs (VELDAT, TEMDAT, SALDAT) 
are included in Appendix 2. Seasonal Profiles of sound speed, 
temperature, and salinity are also included in Appendix 2. 

Representative monthly temperature and salinity profiles 
were requested from HIDAT at Fleet Numerical Weather Central. 
The HIDAT output was used for comparison purposes with output 
from TEMDAT and SALDAT. The comparison was employed to vali- 
date results obtained in this analysis. 

HIDAT salinity and temperature profiles were requested 
for the entire one-degree square, extending from latitude 32°N 
to latitude 33°N and from longitude 79°W to longitude 80°W. 
Typical values for salinity and temperature at standard depths 
for each month were the output for HIDAT. Data for the entire 
square was requested due to the fact that a one-degree square 
is the smallest subdivision possible with the HIDAT system. 


Seasonal averages of salinities and temperatures for standard 
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depths to 100 meters were made of the HIDAT typical profiles 
in order to compare the results with SALDAT and TEMPDAT output. 


The results of the comparison are summarized in Table III. 
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III. RESULTS 


VELDAT analysis of the inshore region (I) of the study 
area produced little in the way of usable results. Data den- 
sity in this region was so sparse as to render the results 
statistically meaningless. Less than 50 stations (all months, 
all years) were observed in region I. 

Results of the analysis of region I exhibited standard 
deviations in the range of 1.2 meters/sec to 8.6 meters per 
second, depending on the season. Information was confined to 


the upper 10 meters of the water column as there were no ob- 
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servations below this depth. Computations were based on a 
maximum of 16 samples for any one season. Table IV is a 
summary of VELDAT results for this region. 

The results indicated for region I clearly show that lack 
of adequate samples of historical data in the form of archived 


STD and Nansen Cast observations precludes use of this method 


to obtain velocity correctors for echo sounding. 

In-situ observations would have to be made of temperature, 
salinity and depth in this region in order to obtain velocity 
corrections. This result is important in itself since it 
identifies a region lacking in data. 

Results of the analysis for region II are exhibited in 
appendix 2. A plot of all stations observed is included as 
appendix 3. Stations are plotted on a seasonal basis. The 
results produced indicated that there was little information 


with which to work for depths greater than 100 meters. 
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| 
| However, the major portion of region II exhibits depths in 
the range of 18 to 95 meters (exceptions along the shelf 
break noted). 

Winter mean sound speeds computed ranged from 1521.5 m/s 
at the surface to 1517.9 m/s at a depth of 100 meters. Standard 
deviations of sound speeds ranged from 4.4 m/s (at 75 meters) 
to 5.7 m/s (at 100 meters). 

Mean temperatures for the same period ranged from 19.24°c 
at the surface to 17.46°c at 100 meters. Standard deviations 
ranged from 1.89°c at the surface to a minimum of 1.53°c at 50 


meters. 


Mean salinities exhibited a spread that ranged from 36.29°/,, 


at the surface to 36.24°/,, at 100 meters with standard devia- 
tions showing a maximum spread of from .10°/,, (surface) to 
.32°/,. (20 meters). 

Results for the spring showed less variability. Mean 
sound speeds for the spring ranged from 1524 m/s at the sur- 
face to 1514.4 m/s at 100 meters. Standard deviations were 
in the range of 4.2 m/s at the surface to 2.6 m/s at 75 meters. 

Mean temperatures for this season ranged from 20.63°C at 
the surface to 16.25°C at 100 meters. Standard deviations 
were in the range of 1.75°c at 10 meters to .87°c at 75 meters. 

Mean salinities ranged from 35.72°/.,,. at the surface to 
36.19°/,,. at 100 meters. Standard deviations of .77°/.,. 
(surface) to .15°/.,, (50 meters) were observed. 

The summer months exhibited the smallest variation in 


sound speeds only in the upper 10 meters. 


Summer mean sound speeds ranged from 1541.7 m/s at the 
surface to 1515.6 m/s at 100 meters. Standard deviations 
ranged from 1.7 m/s (surface) to 11.1 m/s (100 meters). 

Mean temperatures ranged from 27.74°c at the surface to 
16.94°c at 100 meters. Variability as exhibited by standard 
deviations ranged from .78°c (surface) to 3.58°c (100 meters). 

Mean salinities ranged from 35.60°/,, at the surface to 
a maximum of 36.18°/,. (50 meters). Standard deviations 
ranging from .53°/,. (surface) to .21°/,,. (50 meters) were 
noted. 

Autumn mean sound speeds ranged from 1533.9 m/s at the 
surface to 1527.8 m/s at 100 meters. Standard deviations 
ranging from 3.3 m/s (30 meters) to 13.8 m/s (100 meters) 
were computed. 

Mean temperatures for any particular depth during this 
period ranged between 24.73°c (50 meters) and 21.17°c (100 
meters). Standard deviations were in the range of 1.2°c 
(50 meters) to 4.91°c (100 meters). 

Mean salinities ranged from 36.25°/,., at the surface to 
36.36°/,.. at 100 meters. Standard deviations for any depth 
fell in the range of .10°/,, (20 meters) to .34°/,. 100 
meters). 

The results indicate that sound speeds are least variable 
during the spring. This corresponds to the period of when 
temperatures throughout the column are least variable. 

Near surface summer temperatures (1-10 meters) showed 
the least variation for any particular period or depth but 


variation increased substantially below 10 meters depth. 
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Salinities were least variable during the winter and fall, 
but standard deviations of salinity were always less than 
1°/,... If a two sigma variation were examined, only two com- 
puted salinity standard deviations would be in excess of 
1°/.. (surface and 10 meter depth during the spring). 
Examination of the physical properties and seasonal nature 
of waters in this region based on the R/V GILL data as 
presented by Kuroda and Marland (1973) shows that these results 
are as expected. 
Table V is a summary of the Region II standard deviations 
of sound speed, temperature and salinity on a seasonal basis 
abstracted from the results exhibited in appendix 2. 
The relation between temperature and sound speed is al- 
ready well-documented. Sound speeds calculated with Wilson's 
(1960) equation are more dependent on temperature than on 
salinity and pressure effects combined. Therefore it is : 
reasonable to expect the largest variation of sound speed to 
occur when temperature variation is a maximum. A 1°c tempera- 


ture increase will produce a corresponding increase in sound 


speed of 4.5 m/s. Likewise, a 1°/,. Salinity increase produces 


1.3 m/s increase in sound speed and 100 meters depth increase 


will produce a 1.4 m/sec sound speed increase (Ingham, 1975). 
The results obtained, indicate that sound speed variability 


is much more sensitive to temperature variation than to salinity 


variation. In addition, the results indicate that the greatest 


variation in sound speeds does not occur at the surface but 


at some depth below it, corresponding to the depth of greatest 
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temperature change. The depth of greatest temperature change 
varies between 75 and 100 meters depending on the season. 

This boundary layer between the zone responding to climatic 

and mixing effect and a lower zone of constant temperature 
decrease is usually noted as the thermocline depth. However 
the thermocline in this region is not well-defined. The motion 
of this boundary layer in the vertical leads to the large 
variability of temperatures observed at depths between 75 and 
100 meters. 

The comparison of results achieved during this study with 
the averaged HIDAT results revealed that, in general, HIDAT 
temperature averages were 1°c to 2°c higher than TEMDAT 
mean temperatures. HIDAT salinities agree with SALDAT mean 
salinities to within .5°/,. for almost all cases. 

The temperature variation noted is probably due to two 
factors: 

1) The HIDAT routine used considerably more data in 

the form of XBT observations than the data base 
for TEMDAT contained. 

2) The HIDAT data is more heavily weighted toward 
deeper (farther offshore) portions of the area. 
This is also the region of higher temperatures 
since isotherms parallel the coast and temperatures 
increase offshore. The "weighting" is due to the 
fact that XBT observations are not made in the 
shallow (200 meters near shore coastal regions as 


often as they are at depths greater than 200 meters). 
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The significance of the HIDAT comparison lies in the fact 
~ that it was an independent selection scheme used to check re- 


sults obtained in this study. 
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IV. CONCLUSIONS 


The results of the VELDAT analysis for region I were 
less than satisfactory. Low data density in this region pre- 
cludes the determination of sound speed profiles with any 
statistical certainty. Therefore, the use of historical sound 
speed information for correction of echo soundings is pre- 
cluded in this region. 

In-situ measurement of temperature, salinity and depth or 
sound speed directly would have to be made in this region in 
order to correct echo soundings. Extrapolation of this result 
to other near shore regions on the east coast of the U.S. was 
not done. Nevertheless, a cursory examination of the number 
of stations occurring in the near shore region of Marsden 
Square 116, of which the study area is a portion, indicates 
that low data density in these regions will be a limiting 
factor. 

Results of the analysis for the offshore region (II) allow 
several conclusions to be drawn. 

NOS accuracy considerations require that mean sound speed 
be known to within +4 m/sec. Considering the variability in- 
volved in making measurements of parameters from which sound 
speed is computed, it is the author's opinion that this require- 
ment can logically be interpreted to include 95.5% of the i 


observations or a variability of two standard deviations (20). 
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Table 5 shows the value in meters per second of a one- 
sigma variation at each standard depth computed using VELDAT 
for each season. This table also gives the value of a one 
sigma variation for salinity and temperature at each standard 
depth. 

The results summarized in table 5 indicate that two stan- 
dard deviations at any standard depth would exceed the required 
+4 meter per second value for sound speed, with the exception 
of the upper 10 meters during the summer. 

Further, during the winter, all values for one standard 
deviation exceed the +4 m/sec requirement at all standard 
depths. 

Other seasons show values of one-standard deviation which 
exceed the +4 meter per second requirement for some portions 
of the water column. 

The variability of the historical data exceeds the +4 m/sec. 
criteria at the 20 level, hence the question of whether or 
not historical information is adequate to develop echo sounding 
corrections is answered as negative for this region. 

Data generated as a result of this study indicated that | 
changes in present methods for in-situ determination of sound | 
speeds would be possible. | 

The salinity data presented indicates that this is the 
least variable parameter affecting sound speed. Salinity data 
shows that the greatest standard deviation observed was .77°/.., 
and this was atypical. Standard deviations for the most part 


were in the range of .1°/,, to .5°/oo.: 
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The analysis of Wilson's (1960) equation (Bivins, 1976) 
indicates that a natural variability of 3°/,,. in salinity 
could be tolerated if temperature measurements were made with 
sufficient accuracy (+.1°c). 

If one sigma variation is in the range of .1°/,., to .8°/..; 
it is clear that if this were extrapolated to 20 to include 
95.5% of the values, the variability would not exceed 3°/,,. 
Therefore temperature observations alone, and historically 
determined salinities would appear to meet the stated accuracy 
requirements. Such temperature observations could be made 
with the XBT system, instead of the currently used Nansen 
bottles with reversing thermometers and STD probes with in- 
creased efficiency. 

XBT's are presently available with resolution capabilities 
of .01°C in temperature and 0.9 meters in depth (McDowell, 
1978). This increased resolution over previous devices re- 
sults from improved recorder design (McDowell, 1978). 

Accuracy limits of the XBT have been quoted as "less than 
.2°C error 95% of the time” (Naval Oceanographic Office, 1978). 
Accuracy of the depth measurements obtained by the XBT are 
stated as less than 5 meters (Naval Oceanographic Office, 1978). 

It has been suggested that accuracy improvements in XBT 
measurements can be obtained by utilizing an improved "rate 
of fall equation" for the probe (McDowell, 1978). Improved 
factory selection of thermistors used in XBTs including cali- 
bration could also increase repeatability of the measurements. 
Further testing and definition of the accuracy that can be 


obtained utilizing the XBT would be desirable. 
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Finally, the data base used in this study indicated that 
near-shore data for the study area was sparse. This fact may 
be characteristic of NODC files for near-shore coastal regions. 
Further efforts to improve the NODC data base would result in 
an improved determination of variability in the region. 

Investigation of the T-S relation in these regions may 
also prove useful for corrector determination. Where the T-S 
relation proves sufficiently constant, salinity determination 
may be made on the basis of water temperature alone. Further 
study will be necessary to define this relationship and deter- 
mine applicable regions. 

The results obtained during this analysis, while not 
statistically rigorous, present a general picture of sound 
speed variability within a shelf region and answer the question 


of applicability of historical information to the correction 


of echo soundings. 
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Station Location Plots 
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